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Abstract 

Using a density matrix formulation for the effective action, we obtain a set of macroscopic 
equations that describe the spin accumulation in a non-homogeneous ferromagnet. We give a new 
expression for the spin current which extends previous work by taking into account the symmetry of 
the ferromagnetic state through a careful treatment of the exchange term between the conduction 
electrons and the magnetization, i.e., d-electrons. We consider a simple application which has been 
discussed previously and show that in this case spin accumulation is an interface effect confirming 
earlier results arrived at by different methods. 
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Spin-momentum transfer (SMT), was predicted by Berger— and Slonczewski^ and ob- 
served experimentallj*^!^ and studied theoretically by various groups .^i^i'^'i^i^i^^i^'^i^^i^^ . The 
recent work by Zhang, Levy and Fert (ZLF)^ argued that besides the spin torque, the sd- 
exchange term between the conduction electrons and the magnetization gives rise to an extra 
effective field which also can induce switching by precession. They showed that its effect 
can be important to distances larger than l.Onm from the interface. Calculations by Stiles 
and Zangwilii^ and by Slonczewski^ show indirectly that the effect of this field is almost 
absent. These calculations give the impression that a macroscopic treatment, such as the 
ZLF calculation, is either not applicable or simply gives a different result for some unknown 
reasons. However a macroscopic approach is very appealing since it can easily be inte- 
grated within micromagnetics and keeps the problem of spin accumulation within the reach 
of classical methods. In fact as we show elsewhere^, the treatment of spin accumulation 
through a spin density vector provides the most straightforward extension of the Valet-Fert 
theory in CPP structures to non-collinear configurations of the local magnetization.— We 
integrate the Boltzmann equation to obtain macroscopic equations of motion for the spin 
accumulation. In this communication, we show how to extend the ZLF treatment and give 
generalized equations for the spin diffusion in the presence of non-uniform magnetization 
which may be useful for problems involving domain walls. 

As an application of our results, we show that in the uniform magnetization case, the 
effective field predicted by ZLF is indeed vanishingly small for distances larger than l.Onm 
from the F/N interface. This is in agreement with the other calculations mentioned above.-^*^ 
Hence, at least in the uniform case, the macroscopic treatment also predicts that SMT is a 
surface effect. Our equations differ from those derived previously by ZLF by including the 
effect of the magnetization on the electron propagators at the microscopic level. We start 
from a microscopic description of the conduction electrons and the ferromagnetic medium 
and take the semi-classical limit to derive equations for macroscopic quantities in the diffusive 
limit, generalizing previous resulta^*^. We believe that this generalization is necessary for 
metallic elements. Details of the derivation which is based on a path integral approach are 
treated elsewhere.— 

We first introduce the notation: With a spin vector Sj at each lattice point i, the macro- 
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scopic spin vector for the medium is 

N 

S (r) = 5^ (r - r,) . (1) 

1=1 

The interaction between the electrons and the locahzed spins is taken to be of the s-d type, 
of the form 

Hsd = -^jd^ (^) ^* (^)) ■ S (x) (2) 
where J is a couphng constant of the order of 0.1 eV and is a vector whose components 
are the Pauh matrices, 

= 2ze'^Vfc. (3) 
e^^^ is the antisymmetric unit tensor. The Hamiltonian of the theory has the form 

n=Y. /rfrC(r)|^P^(r) + f/(r)-a.BU„(r) (4) 

2 3 

E E (r) ^S*/^ (r) S^ (r) + Hm 

a,P=l i=l 

where B is an external magnetic field, U (r) is a spin-independent potential that includes 
the electric field and TYm is the Hamiltonian of the magnetic system alone without the 
conduction electrons. In the above, we are using units such that h = q/ib = 1- We use 
methods of non-equilibrium field theory to extract the semi-classical equations of motion for 
the spin accumulation m.M^ Our treatment is similar in spirit to that of Brataas, Nazarov 
and Bauer— except that here we use the sd-exchange model to simulate the spin momentum 
exchange between the electrons and the magnetization. Recently Millsl2' and before that 
Berger^ gave a detailed treatment of this exchange in the ballistic regime. Here we focus 
on the diffusive regime which is applicable in large devices such as GPP recording heads 
(see Ref. for a further discussion of this point). To get a macroscopic description of spin 
accumulation, we need to derive equations of motion for the two-point propagators of the 
conduction electrons' field 

^-'(x,y) = (r-^(vl/,(x)vl/+(y))), (5) 

i^,y) = {^s (x)v&+(y)), 
Gn' (x,y) = (T(vI/,(x)vI/+(y))), 
gt^' (x,y) = -(M/+(y)M/, (x)). 
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T is the usual time ordering operator and x is a 4- vector with time and spatial components. 

Using the standard tools of field theory we derive the approximate effective action for 
the conduction electrons and the localized magnetic moments (see Ref. Q and references 
therein for more details). From the effective action we obtain equations of motion for the 
spin current and the magnetic moment. We use the true electron propagator taking into 
account the local magnetization in a self-consistent way. This self-consistent treatment of 
the effect of the magnetization on the electrons is the main goal of this work. We define 



9JlL^(x,y) = -aL,^:'|(x,y) , (6) 

to be the conduction electron spin propagator. The classical polarization of the current m 
is found by averaging over the fast degrees of freedom, in the center of mass coordinates and 
using the quasi-particle approximation.^^ 

The spin current ^, a fourth rank tensor, is defined in the usual way 

where t and x are now coordinates in the center of mass. The spin accumulation vector is 
therefore given by 




To simplify the treatment, we assume that the splitting of the conduction electron bands 
is small with respect to the Fermi energy, i.e., the polarization of the conduction electrons 
is small. We have also used a quasiparticle approximation for the collision term. For slow 
variations in time, we have a modified Fick's law and Ohm's law for the spin accumulation 

(t, x) = -7m^^ (x) - ^'^^dx^m^ (t, x) , (9) 

where 7 is the effective mobility^. This is the major result of this communication which 
improves on previous expressions used for spin accumulation in ferromagnets and is valid 
for a general configuration of the local magnetization. The effective diffusion tensor is 

D'=P(x) = D|| (A-i)'^(x), (10) 
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where D\\ is the diffusion coefficient and 



AM 



with the effective local field given by 



1 



tH, rHy 
1 -tH^ 
1 



(11) 



H = B + JM(x). (12) 

T is the (momentum) relaxation time. The diffusion coefficient in Eq. El is a tensor with 
symmetry about the local axis of the magnetization whereas in ZLF the diffusion coefficient 
is a scalar. This will have important consequences as shown below. In steady state, the 
average magnetic moment m obeys 

V d,, [D^'P (x) d,,mP (x)l = —m'' (x) + J [M (x) x m (x)]^ . (13) 

Next, we apply this result to a case similar to that in ZLF. First, we write explicitly Eq. 
^Jfor a magnetization which is a function only of distance x in the direction of the current, 
M = M(x)z, 



^ (Prrij. (x) ^ (P"m^ ix) 
^ dx"^ dx"^ 
Ply da {x) dm^ (x) ^ / ^ ^-P^yA da (x) duiy (x) _ {x) a {x) uiy (x) ^^^^ 

Da {x) dx dx \ ^ D J dx dx Tsj Tsj 

_ ^ d'^m^ (x) d'^niy (x) 



dx"^ ^ dx"^ 

-{D r,^ly\ da{x)dm^{x) ^ 0% da {x) dniy {x) _ niy {x) a (x) m^; (x) , . 



D J dx dx Da (x) dx dx 



where 



a (x) = tJM (x) , (16) 

and the remaining coefficients are given below. If the magnetization is a function of x, y 
and z, the equations are even more involved. In this case, the transverse components of the 
spin accumulation will not decouple from the components along the magnetizationi^. As 
a consequence, the spin accumulation can be enhanced by choosing appropriate gradients 



of the magnetization along the interfaces. If we assume that the local magnetization is 
uniform and along the z-axis, the equations simplify considerably and we find (instead of 
Eqs. 11 and 13 of Ref.y). 



Dp—m^ (x) + D^y—niy (x) = - JMrriy (x) (17) 
-dJ-^^^ + Dp^ruy (^) = ^ + JMm^^ (x) (18) 

with the diffusion coefficients defined by 

_ _ _ ^11 

1 + {tJM{x)) 

and the off-diagonal terms, which do not appear in the ZLF theory, are 

" "l + (rJM(x))' ^ ^ 

In transition metals, the coefficient \Dxy\ ~ 100 Dp (in the bulk) and D\\ is about three 
orders of magnitude larger than D^y for a transition metal such as nickel. These equations 
are strictly valid for the bulk. Any account of the N/F interface is taken care of by suitable 
boundary condition a^^i^^ . The equations are given in their simplest form to compare to 
Ref. y, which similarly does not take account of the interface. The off-diagonal terms are 
absent in all previous works on spin-momentum transfer. However, the above equations are 



similar to equations derived 



electron gas. In Refs. 



19 and 



jy Hirst — and Kaplan for a circularly polarized itinerant 



20l the off-diagonal terms are attributed to exchange stiffness. 
In our case, the non-diagonal terms D^y are also due to an effective exchange between 
the conduction electrons mediated by the magnetization M along the axis. It is these 
exchange coefficients that sets the scale of spin transfer in the problem when tJM >> 1. 
The equations now reflect the symmetry of the ferromagnetic state as required. The different 
terms presented here will have a signiflcant effect on the Valet-Fert theory when extended to 
treat non-collinear magnetizaton in the different layersi^ and on the result reported in ZLF 
about the effectiveness of the ZLF mechanism, i.e., the transverse fleld. The symmetry is 
easily seen by introducing the complex spin accumulation m, 

m (x) = nix {x) — iniy (x) (22) 



a complex diffusion coefficient Dc, 

D, = Dp + iD^y, (23) 

and a complex relaxation time r^, 

- = — + iJM{x). (24) 



The steady-state transverse spin accumulation obeys, 



(25) 



where 

A,' = T,D,. (26) 
The local z-component of the spin accumulation obeys 



dx^ Kdi ' 



(27) 



where Xsdi is the (longitudinal) spin diffusion length which can be 5-100 nm. The general 
solutions for the complex accumulation are of the form: 

m(a;) = Aexp[-a;/Ac] + 5 exp[a;/Ac]. (28) 

The spin accumulation therefore shows in general an exponential decrease (or increase) with 
some oscillations. In this notation, the spin current is, omitting the electric field, given by 
two components, the perpendicular one which is complex and the longitudinal one, 

(29) 



In Ref. |9l the Landau-Lifshitz equation is 

^ = -7oM X (He// + hMREp) - 7oaM x {}Aref x M) + "Mx^ (31) 

a and h measure the relative strength of the Slonczewski term and the transverse term, 
respectively. In ZLF^, the ratio h/a was found to be almost 2 for their Aj = 2.0 nm, the 
transverse spin accumulation length. To avoid some problems with the boundary conditions 
in the solution of ZLF, we solve the same problem for a semi-infinite ferromagnet using the 



ZLF equations and our equations for the same parameters. We assume, as they did, the 
spin current is initially polarized along the vector u = (0, — sin^^, cos^^) by another layer 
with a magnetization parallel to u. Figure 1 shows the original solution of ZLF for the x 
and y components. In Fig. 2 we show our solution for the same boundary conditions, i.e., 
we assume that the spin current is continuous and the transverse components of the spin 
accumulation decay to zero in the ferromagnet. The spin diffusion length is taken to be the 
same in both cases. Both figures are given for A = 10.0 nm, J = 0.10 eV and Tsf = 10^ — 11) 
s which are typical parameters of a soft permalloy material. With these parameters, the 
spin accumulation in the ZLF case seems to be appreciable for distances up to almost lOnm, 
while in our case it vanishes within l.Onm from the interface. Our result is in agreement 
with other calculationsiS and clearly shows that spin accumulation is a surface effect. This 
can be traced to the fact that the polarization of the conduction electrons in a ferromagnet 
diffuse much more slowly in the direction perpendicular to the magnetization as opposed 
to that along the magnetization due to the s-d exchange in our case. It should also be 
observed that in the limit where spin- flip scattering is important, the oscillations in the spin 
accumulation are "overdamped" . This is in contrast to what the quantum calculations of 
Berger-, Mills-^ and in Ref. i21i give which does not take into account spin relaxation effects. 
These latter oscillations are probably an artifact of the free electron band model as shown 
recently in Refjs^. 

In summary, we have written the equations for the spin accumulation in a non- 
homogeneous ferromagnet in a form which reflects the symmetry of the state and obtained 
a non-trivial self-consistent expression for Fick's law. We find that the predictions of the 
ZLF theory are correct only if we neglect exchange effects. This is clearly not the case 
in a transition metal. We have shown that a correct inclusion of exchange effects greatly 
constrain the spin accumulation to be important only near to the surface. This conclusion 
may not be valid when the magnetization is a function of all three coordinates. 
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FIG. 1: ZLF solution for a half-plane geometry. The x and y components are normalized to one. 




FIG. 2: Our solution for a half-plane geometry for the same parameters as in Fig.l (see text). The 
X and y components are normalized to one. Note the difference in scale along the x-axis (nm) with 
respect to fig.l. 
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